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By G. J. Delio 


SUMMARY 

A typical plant process was studied in which the output increased 
with the input until a definite maximum, was reached. Beyond that maxi- 
mum, a sudden drop with an associated inert zone occurred in the output. 
Turbojet -engine operation in the vicinity of stalled compressor perform- 
ance is a process of this type and was considered in detail. 

The output, as it varied about its maximum value, was subjected to 
harmonic analysis to determine signal possibilities, and the results 
were used in the design of the control. A band-pass filter selected a 
particular harmonic function and used it in a prescribed manner as neg- 
ative feedback to maintain the input that produced maximum average 
output . 

This form of peak holding control, which is one method of obtaining 
information from a peak output for purposes of maintaining that peak, 
was studied for the following criteria: maintenance of maximum peak 

average output, minimum duration beyond peak output, quickness in re- 
sponse to reach peak condition, and type of stability. 

The control investigated maintained an average output less than that 
realizable when the input was at the maximum value. The average output 
was directly related to the perturbation amplitude, the magnitude of drop 
occurring at the peak, and the width of the inert zone in the input. It 
also depended on the center frequency and band width of the band-pass 
filter used to measure harmonic content. 


INTRODUCTION 

There are certain processes in which the plant output increases to 
a definite maximum as the input is increased. A further increase in the 
input is followed by a sudden drop in the output to a value that may be 
uneconomical and perhaps damaging. These processes are found in various 
forms: mechanical, chemical, hydrodynamic, electrical, magnetic, and so 

forth. The need for obtaining the maximum output from these plants pre- 
sents a problem in control, and the control of the described plant pro- 
cess leads to a study of a form of peak holding control. 
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Most controls try to maintain the output either constant or at a 
prescribed change using error data (difference between desired output and 
measured output) supplied continuously in an uninterrupted flow. This 
study is concerned with peak holding control of the described process. 

The control is actuated by error data supplied intermittently and not 
necessarily at discrete moments equally spaced in time. The output is 
sampled periodically, and the average input that causes the output to 
be near its maximum value is maintained. » 

An example of the described plant process is the acceleration of 
gas-turbine engines. Compress or -discharge pressure increases with in- 
creasing fuel flow and reaches a maximum as compressor stall is approached. 

At stall, the pressure suddenly drops. Because of the relation between 
pressure increase and accelerating torque, a control that maintains max- 
imum pressure increase would also maintain maximum acceleration, or 
nearly so. 

A control of this process may be obtained by superimposing a per- 
turbation on the input fuel flow to determine the proximity of stall. 

The control would increase the average fuel flow toward the maximum and 
probe for the maximum pressure at the perturbation frequency until the . 

sudden drop in output pressure is sampled. The drop into the undesirable 
region of operation would determine the desired maximum fuel flow anri 
corresponding maximum pressure ratio across the compressor. 

In order to facilitate the study of peak holding control, a function 
generator was synthesized that simulated the essential features of a 
turbojet-engine compressor, in particular, the pressure drop and hyster- 
esis at stall. Then a control was designed that selected the information 
from the simulator output that characterized the peak output and used it 
in a manner prescribed to cause the simulator output to reach and remain 
near its maximum average value. Both the simulator and control were com- 
posed of units of an electronic differential analyzer. 

The general objective of this report is a study of the character- 
istics of the described form of peak holding control. The control has 
not yet been used in experimental turbojet -engine study. This study is 
concerned with only one method of extracting the maximum amount of in- 
formation from a peak output for the purpose of maintaining that peak. 

The criteria sought in this control process are: maintenance of peak 

average output, minimum duration of time in undesirable operation, quick- 
ness in response to command to reach peak condition, least amount of ex- 
cursion in the undesirable region, and type of stability. * 
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PROCEDURE 

The investigation of control possibilities was performed by first 
determining an idealized compressor function from a study of experimen- 
tal pressure responses in the vicinity of the stall region. Then signal 
possibilities were sought by performing a harmonic analysis on the pres- 
sure output. The results of this harmonic analysis were used as a guide 
in a control design. The control was applied to the simulator, and an 
analysis of control behavior was made. 


Analys is 

The compressor function was determined by a study of pressure re- 
sponse and acceleration resulting from fuel-flow changes at constant- 
speed conditions. When the input fuel flow (with a superimposed pertur- 
bation) is near the value corresponding to a maximum output, the result- 
ing sudden pressure drop causes the output to be rich in harmonics. A 
Fourier analysis of the pressure output in this region was used as a 
guide in the selection of usable information for control synthesis. 

Determination of compressor function . - As shown in figure l(a), 
the compressor-discharge pressure drops suddenly when the critical- 
pressure ratio is reached at stall. The speed shows an instantaneous 
change in acceleration. Both changes indicate stalled turbojet behav- 
ior. Their value at the time of stall represents the peak values to be 
maintained. 

From this type of data, the plots of figures l(b) and (c) were made. 
Both the pressure change and acceleration, at constant -speed conditions, 
increase with fuel flow. The increase becomes less as the fuel flow 
reaches the value sufficient to cause stall. At stall, both the 
compressor- discharge pressure and acceleration decrease. The decreasing 
radius of curvature shown on these figures indicates that harmonic con- 
tent in the output should increase as the peak value is reached. 

If the high-frequency fluctuation during stall is ignored and aver- 
age pressure is substituted (fig. l(a)), if the pressure gain with re- 
spect to fuel flow (at constant speed) is assumed linear and to possess 
the same slope throughout the speed range (fig. 1(b)), if surge behavior 
is neglected, and if arbitrary values are used for an inert zone in the 
fuel-flow input, then the idealized function shown in figure l(d) 
results . 

This function describes an input fuel-flow level w^. that produces 
a linear output pressure increase AP at slope m until the stall 
value AP is reached. At this point there is an instantaneous pressure 
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drop Ap. (Symbols are defined in appendix A.) For the case with an 
inert zone in fuel flow, the output follows the path indicated. Thus 
the combination of the inert zone in fuel flow and the pressure drop 
describes a hysteresis loop. 

The idealized function defined for the cases with and without hys- 
teresis describes a plant process that should serve as a rigid test for 
peak holding control. The constant slope up to the maximum value con- 
tains no information that could warn of the impending sudden drop. The 
sudden drop in the output must be sampled for this control to determine 
the maximum value AP. On the other hand, if the slope decreased with 
increasing input (fig. l(b)), the control could hold at some minimum 
positive slope without the need for sampling the sudden pressure drop. 

The size of the inert zone indicated in figure l(d) is arbitrary. 
Its magnitude has not been definitely established. It could be rela- 
tively small if the input were checked at the instant stall is sensed, 
thus minimizing the amount of excess fuel dumped into the combustors. 
This excess, representing delayed combustion, would maintain the output 
pressure after the controlled input fuel flow had been reduced. If fuel 
flow were checked at the value producing maximum pressure, the only con- 
tribution to hysteresis could be from aerodynamic considerations and 
there would be none from the combustion process . 

Signal possibilities . - If, for mathematical convenience, a test 
signal A sin o^t is superimposed on the input fuel-flow level when the 
input level is well below AWf, the output will contain the fundamental 
frequency in undistorted form. When the average fuel flow is near Aw^, 
the output pressure AP will be distorted, as indicated in figure l(d). 
Assuming that the output f(t) is periodic, a Fourier series of the form 

f(t) = Tp + cos nt + b^ sin nt 

n=l n=l 


may be developed. Now assuming that for small excursions into the un- 
desirable region the output depends on the parallel slope N = constant, 
the Fourier analysis of the output yields 


With hysteresis: 


f(t) . 


m “f - f? ("l - * 2 > 


ma + — (cos cp, 
at i 


n=2 


z x sin n<P / . 

^cos n<P 2 - cos n<p^J — — - ^sin n<P 2 


cos q> ) 
sin n<P-^) 


sin <P + 
— 

cos n<p 


n 


n 


3942 


3942 


NACA RM E56B10 


5 


where 

*1 “ Vl 
^2 = ( D 0 t 2 


stall inception angle 
stall recovery angle 


No hysteresis: 

f(t) = 


mx - ^ (l 


4t 


'1 


^ma - cos <P-^ sin <P + 
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(-l n + l)sin + (.^ . 1)cos siSJl? 


n 


For the case with no hysteresis, the average output pressure for a 
particular input fuel-flow level w f becomes 

s ^ (i - -?) 

But since the maximum output pressure is 

AP = mw f + ma sin ^ 
the loss in the average output pressure becomes 



The amplitudes of the various harmonics of f(t) with no hysteresis 
are plotted in figure 2 as functions of the stall inception angle 
or that instant of the fundamental period when maximum pressure is reached. 
Figure 2(a) sho ws that the ratio of average pressure rise to maximum 
pressure rise AP/AP decreases as the increasing average fuel flow probes 
into the stall region. The maximum amplitudes of the even harmonics de- 
crease with increasing order, and the maximum occurs progressively close 
to the quarter-period time. 

In figure 2(b) the amplitude of the first harmonic starts to de- 
crease at the moment stall is experienced and becomes zero when stall 
duration is for one-sixth of the period. A control based on this sign 
reversal could maintain average fuel flow at a value that would cause 
stall to be experienced during one-sixth of every period with a corre- 
sponding loss in average pressure increase. The amplitudes of the re- 
maining odd harmonics become progressively small with increasing order. 
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The mean-square value of each harmonic (fig. 2(c)) computed from 

o n p p 
®0 ^ + b n 

4 + Z_y 2 
1 

decreases with increasing order, and the maximum value for the first 
harmonic occurs closer to the quarter -period time. When either the odd 
or even harmonic is at its maximum value, stall duration is for n/n 

part of the period. Stall inception is at ^1 - ^ part of the period. 

The effect of summing the harmonics is illustrated in figure 2(d). 
The maximum value increases and stall inception approaches the quarter- 
period time. Thus, a control that used some function of the harmonic 
content of the output as negative feedback could maintain the input fuel- 
flow level near its maximum value. The mean-square values of the har- 
monics bear a ratio to each other, as shown in figure 2(c). A partic- 
ular choice of harmonics (fig. 2(d)) affects the remaining parameters in 
a control; that is, the second harmonic is at a maximum when stall is 
experienced for a quarter cycle of the test frequency, whereas the tenth 
harmonic would need amplification but would be a maximum when stall is 
experienced for only a twentieth of a cycle. The corresponding drop in 
average output pressure, however, would be less. 


Method of Control 

The idealized function in figure l(d) was simulated by an amplifier- 
diode circuit using commercially available components of an electronic 
differential analyzer. The simulator was then used as a plant process 
to be controlled. The control used a band -pass filter that selected 
particular harmonics from the output when stall was experienced and used 
them as negative feedback to reduce the fuel flow. 

Simnl ation . - The nonlinear characteristics of the compressor func- 
tion were simulated with the elements shown in figure 3. The first ele- 
ment is a comparator. It compares the input with a reference limit E^, 
which is analogous to the maximum value AP. When the input reaches 
this value, the comparator output suddenly drops by an amount E^, which 

is analogous to the sudden drop Ap following the maximum value. 

The comparator serves as an input to a limiter (element 2) . This 
element has unit negative slope and limits its output to Eg when its 
input drops below -Eg. The value Eg is analogous to the inert zone 
in the input. When the limiter output is fed back positively, thus 
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closing the loop, the comparator output functions as a bistable multi- 
vibrator of modified form, provided the input sweeps through a suitable 
range. It has a hysteresis loop whose position and area are independent 
functions of the three variable parameters Ep, Eg, and Eg. 

The input to the comparator also serves as an input to a limiter 
with positive unit slope (element 3) . Its output is limited at the val- 
ue of maximum pressure rise E 3 . When this output is summed (element 4) 
with the comparator output, a function generator results that compares 
closely with the idealized function of figure l(d) . The actual function- 
generator output is shown in figure 4, and the detailed description of 
the design is contained in appendix B. The simulator describes an input 
fuel-flow change e^ (volts) that produces a linear output pressure 
change e Q (volts) with a slope m until the stall value E 3 is reached. 
Here the output drops by the amount Ei (pressure drop) to an average 
stall value and remains at this value for further increases in e i . For 
decreasing e i , the output follows the path indicated to describe the 
assumed hysteresis loop. 

Method . - A sinusoidal variation in the input, when the input e^ 
is close to E 3 , produces a pulse of duration t and height equal to 
. If a control could manipulate the input in an effort to minimize 
the pulse duration, it would maintain the plant process mostly in the 
desirable region and would yield an average output near the max imum val- 
ue Ej. The information contained in the pulse must be used as quickly 
and efficiently as possible. 

The frequency spectrum of a pulse is shown in figure 5. Most of 
the signal power is localized at the low-frequency end. However, for 
purposes of control, only the increase in the output is of importance. 

The average output value is not required. The fundamental frequency 
(fig. 2(b)) presents difficulties in its use in control. Also, the con- 
tribution from the frequency range above l/t is negligible. Thus only 
those frequencies between the second harmonic and l/t are used in 
control . 

One possible control is shown schematically in figure 6(a). The 
simplified block diagram is shown in figure 6(b). A band-pass filter 
selects the particular harmonic or group of harmonics from the simu- 
lator output. This output is rectified to produce the absolute value 
|G(s)|, which is then amplified by the amount K. The output of the 
absolute -value device K|G(s)| exists only when the peak value of the 
simulator is reached; otherwise, it is zero. 

A forcing bias C, whose integral produces a ramp velocity that 
drives the input toward the maximum value, is summed with the preceding 
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output. The integral of this sum produces a ramp velocity of C units 
per T seconds during the time existing before the peak value E3 is 
sampled. When the peak value is reached, the generated signal produces an 
average absolute value. This value K|G(s)|, present only when the peak 
value is sampled, subtracts from the forcing bias to reduce the ramp 
velocity. The ramp velocity reduces to zero at the time the amplified 
average absolute value of the filter output is equal and opposite to 
the forcing bias. Under this condition the peak value will be held. 

The integrator output is summed with the perturbation A sin o^t, and 
this sum is used as the input to the simulator; thus the loop is closed. 
The controller output which manipulates the input e^ into the simu- 
lator becomes 


using the constant given in figure 6 (a). 

The use of a perfect ideal filter and absolute -value device was 
assumed with the described theoretical control. The imperfect filter 
does not completely separate the signal from the noise environment. 
However, a minimum-width band-pass filter, used in the control, can con- 
siderably increase the signal-to-noise ratio. The imperfect filter, be- 
cause of its selectivity limitations, will also cause deviation from 
theoretical control behavior. Thus, the entire control depends to a 
great degree upon filter characteristics. 

The filter behavior is described by figure 7(a). The filter is a 
commercially available unit having essentially linear phase shift in the 
pass band together with approximately zero-decibel insertion loss. At 
the minimum band width, the insertion loss is one-half. However, there 
exist a delay time and build-up time that are inversely proportional to 
the band width. They depend mainly on the slope of the phase function 
within the pass band; that is, they increase with selectivity (ref. l) . 

Filter responses to input pulses are shown in figure 7(b). The 
ringing amplitude and time decrease with increasing center frequency 
(fig. 7(b), left) and increases with band width (fig. 7(b), right). The 
amplified absolute values are also shown. 

For the purposes of control only the higher harmonics generated at 
the peak are desired. The direct -current and the fundamental harmonic 
yield only information characterizing the region below the peak value. 

The left side of figure 7(c) illustrates the cut-off characteristics of 
the filter. For a center frequency at the second harmonic and minimum 
band width (the fundamental, is chosen at 60 cps), the filter output was 
considerably attenutated (fig. 7(c), left); at the third harmonic the at- 
tenuation was practically complete. The other record (fig. 7(c), right) 
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yields the absolute value of the filter output; when the input is a sig- 
nal produced by the simulator operating near its peak value under con- 
ditions of hysteresis. 


RESULTS AND DISCUSSION 

The particular control described by figure 6 was "set up" on an 
electronic differential analyzer for study and analysis. Various con- 
trol variables were recorded during stable peak holding action, and 
measurements were made of peak; average output and filter-output values 
at the point of "best stability." The effect on stability of such de- 
sign parameters as band-pass center frequency, band width, and loop gain 
were determined; and general control action is discussed. 


Typical. Recordings of Control Action 

Several analog solutions of the control are presented in figure 8. 

An arbitrary forcing bias -C, which could be some control function such 
as speed error in a gas-turbine engine, was integrated to supply a ramp 
velocity (fig. 8(a), trace 3). A perturbation was superimposed on the 
ramp (trace 5) which served as the input e i to the simulator. The 
simulator output is shown as trace A. When the output reached the max- 
imum value E 3 , wave trains appeared at the filter output G(p) (trace l) . 

The filter was adjusted to minimum band width and a center frequency of 
the third harmonic. The amplified absolute value of the filter output 
appears as trace 2. 

Control action can be described as follows: Figure 8(a) is a re- 

cording of important parameters when the control is operating at the 
best stability condition with moderate ramp velocity. Initially there 
is a ramp velocity determined by the forcing bias and integrating rate. 

The moment maximum output is reached, the action of the control starts 
to reduce the ramp velocity to zero. The main stability criterion is 
the vanishing of the integrand; that is, the forcing bias must be equal 
and opposite to the amplified average absolute value of the filter 
output . 

If the gain of the absolute -value device is further amplified, a 
higher forcing bias must be used to maintain the best stability. This 
is illustrated in figure 8(b). The resulting ramp velocity is increased, 
and the maximum output value is reached in less time. However, the in- 
tegrator output becomes more oscillatory. 

The forcing bias need not be equal and opposite to K|G(p)| . When 
it is less than the average absolute value of the filter output (fig. 

8(c)), the control behaves as though it is too sensitive to the sudden 
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drop at maximum output. The resulting filter output reduces the input 
to the simulator faster than the forcing bias increases it. For the 
particular conditions shown in this recording, the maximum output is 
sampled on alternate cycles only. 

When the forcing bias is too large (fig. 8(d)), the control cannot 
arrest the input when maximum output is reached. The sampling of the 
sudden drop in the output reduces the ramp velocity, but it is not re- 
duced to zero. Eventually, the input drives the output completely into 
the undesirable region of operation. 

In figure 8(e), the sudden drop in the output E-^ has been reduced. 
The generated signal at the output of the absolute -value device also de- 
creases. This result indicated a smaller value of forcing bias to main- 
tain best stability. The ramp velocity is considerably decreased. 

The effect of hysteresis is shown in figure 8(f). At the maximum 
output, the ramp velocity not only reduces to zero, but it is also made 
momentarily negative. The recovery from the sudden drop occurs at a 
latter part of the perturbation period and is equal to an output value 
less than maximum. 

In the hysteresis cases, it is possible to obtain two stable peaks 
(fig. 8(g)). This condition is due to the existence of a double value 
in the rectified filter output. There are two points at which the areas 
per unit time under the rectified filter output are equal; one point 
exists for a short -pulse duration, the other for a longer -pulse 
duration. 

The analog solutions presented in figure 8 were recorded with gal- 
vanometer elements. These elements were not critically damped and had 
dynamic behavior that caused overshooting following sudden drops in the 
simulator output. An oscilloscope used in recording two simultaneous 
functions during alternate periods (fig. 9) indicated that the dynamic 
error was due to the recording equipment. 

The upper photograph on the left side of figure 9 includes alternate 
input and output values of the controlled simulator during peak stable 
operation for no hysteresis. The lower photograph is for a condition of 
hysteresis. The pressure drop and recovery from the drop are instanta- 
neous. The duration in undesirable operation (pulse width) is less for 
the case of no hysteresis. 

The upper photograph on the right side of figure 9 illustrates the 
relation between the controlled output of the simulator when it is near 
maximum and the output of the absolute -value device. The lower photo- 
graph is a record of the same functions when the controlled output is 
maintained at the second stable point. The pulse width covers nearly 
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the entire cycle during which the simulator is operating in the unde- 
sirable region. The ringing time of the filter for this case is less 
than the pulse duration; that is, the transient caused by the sudden 
output drop had time to completely attenuate before recovery was caused. 


Measurements of Average Values at Best Stability 

The average simulator output and average absolute -value -device out- 
put were measured at the peak values under conditions of best stability 
with appropriate averaging meters . These values were measured to deter- 
mine the effect of control design parameters on stability. The param- 
eters varied were band-pass center frequency, band width, magnitude of 
output drop at maximum output, width of inert zone in input, perturba- 
tion amplitude, and gain of absolute -value device. 

The data for figures 10 to 12 were obtained when the control was 
acting in its most stable manner. This occurred when the forcing bias 
was approximately equal and opposite to the average output value of the 
absolute-value device. The effect of varying the band-pass center fre- 
quency and magnitude of drop at maximum output E-^ is shown in fig- 
ure 10(a) . The peak average output increases with band-pass center fre- 
quency but decreases with increasing output drop. The output of the 
absolute -value device decreases with increasing center frequency but in- 
creases with increasing output drop. The maximum usable forcing bias 
varies with the output of the absolute -value device. Thus, larger out- 
put drops permit a faster ramp velocity but cause the peak average out- 
put to reduce. 

The effect of varying the filter pass band width is illustrated in 
figure 10(b). The peak average output decreases with increasing band 
width and with the size of the inert zone in the input. The output drop 
is kept constant. The average output of the absolute -value device in- 
creases with increasing band width and inert zone. This permits faster 
ramp velocities to be associated with the wider band widths. 

Hysteresis effects are shown in figures li(a) and (b). The peak 
average output decreases with increasing inert zone and with increasing 
output drop. The average output of the absolute -value device remains 
essentially constant with variation of inert zone but again increases 
with output drop. 

As the perturbation amplitude increases, the peak average output 
decreases (fig. 12). However, the average output of the absolute-value 
device, and hence the ramp velocity for best stability, remains con- 
stant. The time delay and lag of the actual filter, coupled with the 
voltage drop across the diodes contained in the simulator (see fig. 17), 
did not permit the pulse duration to reach zero. Therefore, the pertur- 
bation amplitude must be at least equal to the inert -zone width for the 
control to operate properly. 
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The peak average output of the controlled simulator did not vary 
with the gain of the absolute -value device (fig. 13). The output of 
the absolute-value device varies linearly with its gain and allows the 
ramp velocity to vary accordingly. 


Analysis of Control Behavior 

The control is an example of a servo (including computor) operating 
with pulse data (ref. 2). The control begins operating with an initial 
error (forcing bias). Immediately after a drop occurs in the output, a 
measurement is made and the error is decreased. After the next measure- 
ment, the error is again decreased. Instability can result from either 
overcorrection or undercorrection. Initially, the forcing bias is acted 
upon by only the forward part of the loop. When the peak is approached, 
the discontinuities present in the simulator produce pulse inputs to the 
filter, and the loop is closed. Thus the control depends heavily upon 
the performance of filters with pulse input data. 

•The transfer function of the actual filter used is given in figure 
7(a). If, for simplicity, the filter is assumed critically damped (i.e., 
sin (SQt * ooQt), its indiciad response would be as shown in figure 14 
(ref. 3). This theoretical solution agrees with the measured solution 
recorded in figure 7(b) as do also the absolute values. The dominant 
feature of the response (absolute value) is an oscillation whose period 
is short compared with the important time constant of the control sys- 
tem, the integrating rate. The integral of the absolute value of the 
filter response is given in figure 15. For at < 1, the integral is 
approximately zero; for at > 20, the contribution to the integral is 
approximately zero. Thus the corrective value existing as feedback at 
the time at can be approximated as dead time plus first-order lag. 

This corrective action exists only after a sudden drop in the output is 
experienced . 

The continuous operation of the control loop can be used in study- 
ing stability during corrective action. Initially, the forcing bias 
provides a ramp input that drives the output toward its maximum value . 
Superimposed on the ramp input is the perturbation signal A sin o^t . 

The instant the input reaches a value causing maximum output, the sudden 
drop in the output introduces a step into the filter. This step pro- 
duces the response given by figure 14, whose integral (fig. 15) becomes 
the corrective value. 

For pulses whose duration is short compared with the perturbation 
period, the shape of the absolute value of the filter response resembles 
that due to a step input. Therefore, for simplicity, a step input to 
the filter will be considered. 
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The ramp velocity (fig. 16) initially is 

at 



C d(at) 


and the input to the simulator at any time is the sum of this and the 
instantaneous value of the sinusoidal variation. When this input is at 
the maximum value, a step is introduced into the filter. The integral 
of the absolute value of the filter output is 


This integral acts to reduce the ramp velocity to zero (end of region l) 
At this time, the integrals are equal but opposite. Because the correc- 
tive action exists only during part of the perturbation period, the ramp 
velocity must be made negative (region 2) to ensure that the input re- 
mains below the value corresponding to maximum output. 

The corrective integral (fig. 15) tends to a maximum value; the 
ramp integral is constantly increasing. Therefore, the two integrals 
will again produce a zero ramp velocity (end of region 2) . Thus, zero 
ramp velocity occurs twice during each perturbation period if the forc- 
ing bias is made equal and opposite to the average absolute value of the 
filter output. This is the condition for best stability. 

If the corrective action is insufficient, the control is unstable 
in that it drives the output into the undesirable region (beyond peak 
output, fig. 8(d)). If the corrective action is too great, the control 
will cause operation below the peak output. Peak output will not be 
reached during every cycle of the perturbation period (fig. 8(c)). 

When the corrective action is exactly sufficient, the peak value will 
be reached every perturbation period (fig. 8(a)). Proper operation thus 
constitutes a sampling type of peak holding control. 


A typical plant process was studied in which the output increases 
with input up to a certain maximum value, and further increases cause 
the output to drop suddenly into an undesirable region of operation 
(turbojet -engine operation in the vicinity of stalled compressor per- 
formance). The information was incorporated into a specially designed 
function generator that simulated the process . 



SUMMARY OF RESULTS 
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The simulator was used to investigate the possibilities of a peak 
holding type of control. Sinusoidal variations in the input to the sim- 
ulator resulted in a distorted variation in the output when the output 
was near its peak value. The generated information was used as negative 
feedback to maintain the peak value. 

The control investigated maintained an average output less than 
that realizable when the input was at the maximum value. The average 
output was directly related to the perturbation amplitude, the magni- 
tude of drop occurring at the peak, and the width of the inert zone in 
the input. It also depended on the center frequency and band width of 
the band -pass filter used to measure harmonic content. The ramp veloc- 
ity depended on the feedback path gain and magnitude of drop at the 
peak value . 

The control was stable for all ramp velocities below some optimal 
velocity that depended on both loop gain and magnitude of drop. At the 
optimal setting, the ramp velocity was maximum and the input had minimum 
variation . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 2, T955 


3942 


3942 


* 


MCA KM E56B10 


15 


APPENDIX A 
SYMBOLS 


A 


a,b 

C 

D 1> D 2'- D 3> 1) 4 

E 

E 1 

E 0 


% 


e . 

1 


o 

f,G 

K 


K 

2 

!£' 


G(s) 

,-l 


m 

N 

N 

n 

AP 

AP 


perturbation amplitude, v 

Fourier coefficients 

forcing bias, v 

diodes 

voltage 

voltage simulation of Ap 

voltage simulation of inert zone in input 

voltage simulation of AP 

corrective voltage, v 

feedback voltage, v 

instantaneous input voltage, v 

instantaneous output voltage, v 

functional relations 

feedback path gain 

average absolute value of filter output, v 
inverse Laplace transform 
slope, AP/Awjp 

engine speed, rpm or percent rated 
acceleration, rpm/sec 
Fourier series harmonics 

pressure change from equilibrium, lb/sq in. 
maximum pressure increase, lb/sq in. 
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r 


2T 

Ap 

s 

T 

t 

w f 

AWf 

w f 


M- 

T 


t 1' t 2 

*L 

(p 2 

<°c 

CD 0 


average output pressure , lb/sq in. 
pressure drop at maximum, lb/sq in. 

Laplace operator 
perturbation period, sec 
time, sec 
fuel flow, lb/hr 

fuel -flow change from equilibrium, lb/hr or percent rated 

input fuel-flow level 

fuel flow corresponding to AP 

filter parameter, radians/sec 

amplifier gain 

integrator time constant 

0 ) c / 2 ji, sec 

band-pass -filter variables 
o^t, radians or deg 

stall inception angle, radians or deg 
stall recovery angle, radians or deg 
band -pass -filter center frequency, radians/sec 
perturbation frequency, radians/sec 
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APPENDIX B 


COMPRESSOR SIMULATION 

A typical engine was studied for its pressure response to fuel- 
flow changes with particular emphasis on its stall behavior. The sharp 
pressure drop at stall (fig. 1 (a)), the associated dead zone during stall, 
and the hysteresis loop that the combination implies was duplicated by 
an amplifier -diode circuit. A comparator coupled to a series output 
limiter was connected in a closed loop to form a modified bistable mul- 
tivibrator, and the response of the multivibrator was summed with an- 
other series output limiter. The result was a function generator that 
approximated, in a simplified manner, the dynamic behavior of the 
compressor. 

The comparator reference voltage duplicates the stall point, and 
the limiter determines the size of the dead zone. The comparator also 
determines the amount of the pressure drop. These nonlinear character- 
istics were reproduced by amplifier-diode circuits. The circuits used 
diodes having a high back impedance to produce abrupt changes in slope 
that depend on the magnitude of one or more voltages . Diodes rather 
than relays were used because they are faster and usually require less 
equipment. Diodes also permit the use of stabilized direct -current 
computing amplifiers. The error introduced by diode use is given in 
figure 17. 

The analyzer circuit that simulates the instantaneous pressure 
drop at stall is the comparator (fig. 18(a)). A variable voltage Ej_, 

which is a measure of the pressure drop at stall, figure 1 , can be 
turned on or off at a coincidence of two other variable voltages. Eg 
representing the pressure increase at stall and e^ representing the 
input fuel flow. 

The comparator is a switch that operates when the input voltage 
e^ equals the negative of the reference voltage (e^ - E 3 ) . If e^ = 0, 

increasing positively, diode D 2 conducts and clamps the output to 
zero (within the error specified in fig. 17). When e^ = E 3 - e^, a 
transition is initiated and D 2 stops conducting. Then, for a further 
incremental increase Ae^ = E, the full amplifier gain produces an out- 
put e Q = -| 1 E. The output increases negatively until -|xE = -E^. At 
this value D^ conducts and clamps the output to -E^ . Further in- 
creases in the input produce no further change in the output. 
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If = 0, increasing negatively, Dg conducts and clamps the out- 
put to zero. Thus the comparator has zero gain in either of its two 
stable states. The output must be either zero or -E^ according to 
whether the sum of the input and reference is negative or positive. 
During the transition, the output slope is p (70,000 for the amplifiers 
used) . 

In order to simulate dead zone (inert zone in fuel flow), a series 
output limiter is coupled to the comparator output and the limiter out- 
put is fed back to the comparator input to form a closed loop (fig. 18). 
The inert zone in fuel flow (approximated from such data as fig. l) is 
then determined by Eg. 

A series output limiter is defined as an amplifier -diode circuit 
employing one or more diodes in the feedback path that decrease the 
current through the summing node by diode conduction. It produces nor- 
mal scale changing until some limiting value of the output is reached. 
For the given circuit (fig. 18(b)), the output slope is negative. If 
the input is positive, decreasing, the output is negative. As the input 
passes through zero, increasing negatively, the output passes through 
zero, increasing positively. When the input reaches a value for which 
the output is Eg, Dj conducts and thereby clamps the output to Eg 

for increasingly negative inputs . 

The coupling of the series output limiter and comparator to form 
a closed loop (fig. 18(c)) forms a modified version of the well-known 
bistable multivibrator. In this circuit there is a hysteresis loop 
whose position and area can be independently adjusted; that is, the 
stall pressure Ej, the pressure drop E^, and the inert zone Eg can 
be independently adjusted to the values determined from engine data. 

If the input voltage is initially negative, increasing positively, 
and the feedback voltage e^ is momentarily zero, Dg conducts to clamp 

the output to zero. The output remains zero until the input is equal 

and opposite to -Ej (e i - Ej = 0). At this value the transition to the 

other stable state starts. Then Dg stops conducting. The output var- 
ies scalarly, with slope -p, with the input until it reaches the value 


The series output limiter provides positive feedback e^ to the 
comparator. During the transition the comparator output is limited to 
the value Eg. The feedback voltage does not disturb the comparator 

because it is positive and thus does not cause the value at the summing 
junction of the comparator to go through zero. Therefore, the output of 
the comparator remains at -E-|_ for further increases in the input. 
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As the input reverses direction (becoming less positive), the out- 
put remains clamped to -E^ until e^ + e-^ - = 0. The comparator 

thus switches at the lower value E 3 - to zero output. The output 
remains clamped at zero for all values less than E 3 - e^. After the 

transition, the feedback voltage again is zero, thus allowing the com- 
parator to switch, in increasingly positive inputs, at the original val- 
ue Ej. Note that for all | E-^ | < | I the hysteresis loop remains 
oa square because the feedback voltage will not be limited by D 3 for this 

case. 

to 

In order to simulate pressure response to fuel-flow changes, a 
slope (fig. 19(a)) must be added to the multivibrator. As shown in the 
figure, the output slope is positive. If the input to the series out- 
put limiter is negative, the output is negative, and both will be zero 
at the same time. For an input increasingly positive, the output will 
also increase positively until it reaches a value E 3 . At this value 
conducts and clamps the output at E^ for further increases in the 
input . 

When the outputs of the limiter and comparator are summed, a func- 
tion generator is formed that simulates compressor dynamic behavior. If 
the comparator switches at E 3 and the limiter also limits at Ej, then 
for any set values of Ej_, Eg, and E 3 an independent variable serving 
as a simultaneous input to both the comparator and limiter results in an 
output function as shown in figure 19(a). Specifically, if the input is 
fuel flow (lOO v = 10,000 lb/hr), the output is the negative of the 
compressor -discharge pressure (100 v = 10 lb/sq in.). The slope at 0.01 
(lb/sq in. )/(lb/hr) would be the pressure fuel-flow gain at constant 
speed (fig. 4=) . 

The theoretical, function-generator performance is shown graphically 
in figure 19(b). The proposed function generator was set up on a dif- 
ferential analyzer, and the actual output of the simulator is shown in 
figure 3. The voltage drops across the diodes, as shown in figure 17, 
are in evidence. Under the worst conditions the departure from the 
ideal response is on the order of 1 volt. When referred to a machine 
unit of 100 volts, this represents an error of 1 percent. 
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APPENDIX C 


INSTRUMENTATION 

Three separate recorders were used for solution recording: 

(1) A continuously recording oscillograph to record time histories 

(2) An oscilloscope to record instantaneous functions at two suc- 
cessive periods 

(3) A function plotter to record the steady-state relation between 
an input and output function 


Oscillograph 

A standard oscillograph with galvanometer elements and light- 
sensitive paper was used to record time histories. The computing ampli- 
fier rating of 5 milliamperes and range of 100 volts dictated a minimum 
load impedance of 20,000 ohms. Galvanometer elements with a damped sen- 
sitivity of 2 inches per milliampere were used to provide a gradient of 
40 volts per inch on the paper. 


Oscilloscope 

A standard, 5-inch, single-beam oscilloscope was used with a blue 
fluorescent tube having minimum phosphorescence. In conjunction with a 
Polaroid camera, transients were recorded that had a duration of one 
period of the disturbance function. Two successive periods were dis- 
played on the single-beam cathode-ray tube by using an electronic switch 
and adder. 

A sweep generator (fig. 20 ) operating at the test frequency of the 
control triggered a flip-flop circuit. The Eccles-Jordan trigger cir- 
cuit driven from the sweep generator had cathode -follower outputs and 
produced a 30-cps square wave, that is, half the basic repetition rate. 

Two particular solutions, occurring at successive periods (figs. 7 
and 9), were displayed on the scope by the arrangement shown. The sig- 
nal input from the switch turned on each of a pair of diodes on alter- 
nate cycles . Each of the two inputs was connected to a diode through a 
unit -gain (variable) amplifier. A signal from the other side of the 
switch was added to the two diode outputs. An adjustable resistor (sep- 
arator) caused the displays from the two input signals to separate ver- 
tically on the cathode-ray tube. The signal path through the switch was 
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direct -coupled , thus the signals were switched provided the separating 
resistor was adjusted for zero separation. However, there was some 
noticeable switching transient. 


Function Plotter 

The function plotter was commercial equipment that had a writing 
speed of 12 inches per second on both axes. This plotter produced such 
figures as figure 4. 
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(a) Response of turbojet engine to sudden increase in fuel flow. (Calibration: 
compressor -discharge pressure, 6.96(lb/sq in.)/in.; turbine -discharge pres- 
sure, 8.3(lb/sq in. )/in. ; engine speed, 200 rpm/in.; small-slot fuel pres- 
sure, 17l(lb/sq in.) /in.) 

Figure 1. - Typical turbojet -engine performance. 
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(b) Pressure change due to fuel-flow change at constant speed 
in vicinity of stall. 


Figure 1. - Continued. Typical turbojet-engine performance. 
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(c) Variation of acceleration with fuel flow. Engine speed, 70 percent rated. 

Figure 1. - Continued. Typical turbojet -engine performance. 
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No hysteresis 
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Hysteresis 


f(t) 
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f(t) = raw f + mA sin oc^t 
(0 < t < t x , t 2 < t < T) 

f(t) = mw f + mA sin a^t - Ap 

(t^ < t < t^) 


(d) Idealized compressor response to forcing fuel flow, 
speed. 


Constant engine 


Figure 1. - Concluded, Typical turbo jet -engine performance. 
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Figure 2. - Harmonic content of function f(t) with no 
hysteresis . 
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Figure 2. - Continued. Harmonic content of function 
f(t) with no hysteresis. 
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(c) Mean-square value, ^Stall duration is rt/n for maximum a n or b n , 
at zero ramp velocity, and stall inception is at it/2 ^1 - 

Figure 2. - Continued. Harmonic content of function f(t) with no 
hysteresis . 
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Figure 2. - Concluded. Harmonic content of function 
f(t) with no hysteresis. 
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Figure 3. - Block diagram of simulator. 
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Figure 4. - Actual output of simulator. 


E 3 = 50.0 v; E 2 = 20.0 v; E 1 = 30.0 v. 


CM 


NACA RM E56B10 


Relative signal power 


32 


NACA BM E56B10 



Frequency 


Figure 5. - Frequency spectrum of pulse. 
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e i = 20/(K|G(p)| - C )dt - A sin CDQt 


(a) Problem board wiring. 

Figure 6. - Peak-holding-control loop. 
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Filter and absolute-value device 
(b) Simplified block diagram. 

Figure 6. - Concluded. Peak-holding-control loop. 
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(a.) Filter frequency response. 

Figure 7. - Band-pass filter characteristics. 
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(b) Filter responses to input pulses. 
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(5) Absolute value of filter output 

(6) Filter input for (5) (420-420 cps) 

(7) Absolute value of filter output 

(8) Filter input for (7) (560-560 cps) 


(c) Cut-off characteristics of filter, left; filter responses 
at peak conditions, with hysteresis, right. 


Figure 7. - Concluded. Band- pass filter characteristics. 
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Figure 8. - Continued. Analog solutions of peak holding control. 
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Figure 8. - Continued. Analog solutions of peak holding control. 
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(f) Effect on stability due to hysteresis. E-j_ =* 30 v. Eg , 25 v; E3, 50 v. 
Figure 8. - Continued. Analog solutions of peak holding control. 
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(g) Existence of two stable points. E^, 30 v; E 2 , 10 v; E 3 , 50 v. 
Figure 8. - Concluded. Analog solutions of peak holding control. 
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(2) Simulator output, no hysteresis 

(3) Simulator input 

(4) Simulator output, hysteresis 



(5) Absolute value, filter output 

(6) Simulator output, no hysteresis 

( 7 ) Absolute value, filter output 

(8) Simulator output, 2 n ^ stable point 


Figure 9. - Simulator and filter outputs at stable conditions. 
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(h) Variation of filter band width. 
Pressure drop at stall, E]_ = 50 v, 


Center frequency o> c = IOodq, 


Figure 10. - Concluded. - Effect of filter parameters on measured 
system variables. Fundamental frequency, cjoq = 120it; pressure at 

stall, Ej = 50 v; sinusoidal amplitude, A = 10-y/i2 v; gain in feed- 
back path, K = 20. 
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(a) Average output pressure. 


Figure 11. - Effect of hysteresis on measured 

E 3 = 50 v; sinusoidal amplitude, A = 10-/2 
minimum hand width. 


system variables. Fundamental frequency, = 120*; pressure at stall 
v > 6 & in in feedback path, K = 20; filter center frequency, cd c = 10u)q; 
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Figure 11. - Concluded. Effect of hysteresis on measured system variables. Fundamental frequency, cjdq * 120jt; pres- 
sure at stall, E 3 = 50 v; sinusoidal amplitude, A = 10-^2 v; gain in feedback path, K = 20; filter center frequency, 
o> c = IOcjoq; minimum band width. 
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figure 12. - Effect of sinusoidal amplitude on measured system variables. Fundamental fre- 
quency, coq = 120*; pressure drop at stall, E x = 50 v; hysteresis, E 2 * 0; pressure at stall, 

Ej 3 50 v; gain in feedback path, K = 5; filter center frequency, cu c = IOuoq; filter band 
width, 2 ocq to 2Oo>0. 
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Figure 13. - Effect of feedback path gain on measured system variables. Fundamental frequency, odq = 120jt; hysteresis, 
E2 = 0; pressure at stall, E3 = 50 v; sinusoidal amplitude, A = 10-^2 v; filter center frequency, co c « IOcdq; minimum 
filter band width. 


2^62 


NACA RM E56B10 


Relative amplitude 


CE-7 back 


3942 



Figure 14. - Indicial response of filter. 
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Figure 15. - Time integral of absolute value of filter indicial response. 
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Figure 16. - Theoretical behavior of peak holding control. 
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( a) Comparator. 
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(b) Series output limiter. 



Figure 18. - Modified bistable multivibrator possessing a hysteresis inn-n ... 

adjusted. R, 0.1 megohm. ^ loop whose P° sition and area can be independently 


Cn 

CD 


MCA RM E56B10 


01 

CD 





Figure 19. - Theoretical function generator that simulates compressor behavior at stall at constant-speed 
conditions . 
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Figure 20. - Arrangement for independent gain adjustment and simultaneous display of 
two solutions . 
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